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We have measured the current-voltage characteristics of a single Josephson junction placed in 
a high impedance environment. The transfer of Cooper pairs through the junction is governed 
by overdamped quasicharge dynamics, leading to Coulomb blockade and Bloch oscillations. Exact 
duality exists to the standard overdamped phase dynamics of a Josephson junction, resulting in 
a dual shape of the current-voltage characteristic, with current and voltage changing roles. We 
demonstrate this duality with experiments which allow for a quantitative comparison with a theory 
that includes the effect of fluctuations due to finite temperature of the electromagnetic environment. 

PACS numbers: 73.23.Hk, 73.40.Gk, 74.50.+r 



Duality often plays a central role in our understanding 
of the physical world. Some physical models, for example 
the harmonic oscillator, contain the remarkable property 
of self duality, where we can map the model back onto it- 
self in such a way that the role of physically conjugate de- 
grees of freedom is interchanged. Another generic model 
which exhibits self duality is that of a quantum particle 
in a cosine potential coupled to an ohmic heat-bath lj. 
This model also describes the electrodynamics of a cir- 
cuit consisting of a single Josephson junction shunted by 
a resistor |2J. In this Letter we present an experimental 
verification of this self duality. 

In the standard theory of the Josephson effect 0, Q , a 
current / which is less than the critical current Ic flows 
with zero voltage difference across the tunnel junction, 
corresponding to a static state where the phase differ- 
ence across the junction has no time dependence. This 
region of the current- voltage characteristic, known as su- 
percurrent branch, is well-separated from the dissipative 
quasiparticle branch by a voltage of twice the supercon- 
ducting gap V — 2A/e. However, since early work by 
Ivanchenko and Zil'berman p|, it is known that in the 
presence of an impedance Z env {uj) close to a Josephson 
junction with small capacitance C, phase diffusion will 
occur, and therefore dissipative behavior exists also for 
the supercurrent with finite, time averaged voltages (V) 
below the gap. The resulting DC current-voltage charac- 
teristic (IVC) of the junction consists of a supercurrent 
peak, with a tail at finite voltages originating from the 
Josephson oscillations at the frequency fj=2e{V)/h. 
The detailed IVC of the tunnel junction depends on the 
charging energy Eq = e 2 /2C, the Josephson coupling 
energy Ej — HIc/2e, and the impedance of the envi- 
ronment Z emr {uj). When the impedance is frequency in- 
dependent 5fte[Z cnv ] = R, phase diffusion occurs if Q = 
tt(R/Rq)^Ej/2E c < 1, where R Q = h/4e 2 = 6.45 kQ 
is the quantum resistance. This limit, usually referred 



to as overdamped phase dynamics, was analyzed in |5|, 
where the IVC of the junction was calculated in the pres- 
ence of thermal fluctuations due to Nyquist noise in the 
resistor. Although this analytic result had been known 
for quite some time, the supercurrent peak was only re- 
cently measured for a small-capacitance Josephson junc- 
tion embedded in a carefully designed low impedance en- 
vironment, where R <C Rq for all frequencies g|. 

The picture developed in Q and experimentally veri- 
fied in |6j is based on the classical, non-linear phase dy- 
namics of the junction and is valid in a low impedance 
environment. When R ~ Rq quantum fluctuations of 
the phase become important and a theoretical descrip- 
tion of incoherent tunneling of individual Cooper pairs 
can be applied to calculate the junction IVC for an arbi- 
trary Z cnv when Ej <C Eq ■ This theory predicts that 
the supercurrent peak moves to higher voltages, opening 
a Coulomb gap in the IVC of the junction. 

Of great interest however is the case when Ej > Ec 
and R » Rq, usually referred to as the underdamped 
case (Q 3> 1), where the incoherent tunneling picture 
breaks down. Here a situation exactly dual to the above- 
mentioned classical overdamped phase dynamics occurs. 
Based on duality arguments [a and a quantitative the- 
ory in terms of Bloch bands [SjllGj, the IVC of the junc- 
tion was predicted to show a voltage peak near zero cur- 
rent corresponding to Coulomb blockade, and a tail at 
finite currents corresponding to Bloch oscillations with 
frequency Jb — (I)/2e. During each period of the os- 
cillations exactly one Cooper pair tunnels through the 
junction. In this picture the junction circuit is described 
by the classical, non-linear dynamics of the quasicharge, 
which is the conjugate variable to the Josephson phase. 
In view of duality, when the phase dynamics is under- 
damped the corresponding quasicharge dynamics is over- 
damped. The quasicharge then evolves in time according 
to the Langevin equation dq/dt = {lb + 51) —V/ R. Here, 
lb is the bias current and 51 a random noise component 
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FIG. 1: (A) and (B): SEM images of the SQUID and single 
junction sample configuration. A symmetric bias is applied 
from the top-left array to the bottom-right array, (a) and 
(b): IVC of the current bias lines for the two configurations, 
(a) from top to bottom: two SQUID arrays in series in the 
superconducting (Rq » 50 kS7) and insulating (Ro ~ 50 MJ7) 
case, (b): two non-SQUID arrays in series (Ro » 10 MJ7). 



induced by the resistance; the voltage across the junction 
V = dEo(q)/dq is given by the derivative of the lowest 
Bloch energy band (see inset of Fig. 2). The IVC of a 
single junction can be analytically calculated from the 
Langevin equation for the quasicharge [9|, llOl lllj . 

Few experiments have probed a single Josephson junc- 
tion in the regime R 3> Rq due to difficulties in de- 
signing an environment with such high impedance at 
the frequency fs (~ 10 10 Hz). Small on-chip resistors 
[l2l Il3| and two-dimensional electron gas beneath the 
junction [l4| have been used to achieve a high impedance 
environment. In our experiment, Josephson junction 
SQUID arrays are used to bias the junction. The ad- 
vantage of the SQUID configuration is that the effective 
impedance of the environment can be tuned in situ by 
applying a magnetic field perpendicular to the SQUID 
loops. Although no systematic characterization of the ar- 
rays impedance at frequencies of the order of fs has been 
performed, they have been successfully employed |l5j 
to demonstrate how the environment induces Coulomb 
blockade and Bloch oscillations in a single Josephson 
junction in the weak coupling regime Ej < Ec- 

In this paper we present experimental results on sin- 
gle junctions where 5fte[Z env ] 3> Rq. We have studied a 
single junction with SQUID geometry (Fig. 1A), which 
allows for a systematic study of the IVC of the same sin- 
gle junction as it is tuned from strong (Ej > Ec) to 
weak coupling (Ej < Ec), with the magnetic field. We 
also studied a single non-tunable junction (Fig. IB) with 
strong coupling, where the exact dual of the overdamped 
Josephson effect is realized. Here we make for the first 
time a detailed quantitative comparison with theory. 



The Al/AlOx/Al tunnel junctions are fabricated by 
double angle evaporation through a mask patterned by 
electron beam lithography. The samples are mounted in 
a RF-tight copper box and measured in a dilution refrig- 
erator with base temperature of 15 mK. No special cold 
microwave filters were implemented in the cryostat, as 
the arrays themselves are acting as extremely good fil- 
ters, protecting the CPT from fluctuations generated in 
the bias circuit. The single junctions are biased by four 
Josephson junction arrays and the IVC is measured in a 
four point configuration. One pair of arrays is used to 
apply a symmetric bias; the current is measured directly 
with a current preamplifier (modified Stanford Research 
System SR570). The voltage across the single junction is 
measured through the other pair of arrays with a high in- 
put impedance differential amplifier having an input cur- 
rent of 3 fA (Burr-Brown INA116). Since the impedance 
of the voltage leads can be in the Gfi range, even such 
small input current can lead to voltages comparable with 
the blockade voltage of the single junction, making the 
measurement of the voltage over the junction impossible. 

In the sample layout presented in Fig. 1A, the single 
junction has a SQUID geometry, and it is biased by four 
SQUID arrays. Each array consists of 60 SQUIDs with 
a nominal area of 0.06 /1m 2 . Each junction of the cen- 
tral SQUID has an area of « 0.02 /zm 2 . The loop of the 
single junction SQUID is designed to be 10 times larger 
than the loop area of the SQUID arrays, enabling a pe- 
riodic modulation of the Ej/Eq of the single SQUID 
which is incommensurate with the period of modulation 
of the environment impedance. As the magnetic field is 
increased, the SQUID arrays undergo a superconducting- 
insulator quantum phase transition |l6j . This results in 
an increase of the zero bias resistance Rq of the arrays 
over several orders of magnitude (from kil up to GQ). 
Figure la shows the IVC of the two biasing SQUID ar- 
rays at magnetic fields corresponding to the maximum 
and minimum value of Rq. 

Figure 2 shows the IVC of the single junction SQUID 
at two different values of magnetic field. The field values 
were chosen so that the zero bias resistance of the bias- 
ing leads was the same for both curves (Ro w 10 Mfl). 
However, the Ej /Ec ratio of the single junction SQUID 
is different for the two curves. Curve A in Fig. 2 corre- 
sponds to the maximum value of the Josephson coupling, 
while in curve B the Josephson coupling was at a min- 
imum. We estimate the values of Ej/Ec for the two 
curves from the samples parameters. The charging en- 
ergy Ec = e 2 /2C = 45 [ieV is calculated from the junc- 
tion area, which gives a capacitance C ~ 1.8 fF assum- 
ing a specific capacitance of 45 fF //1m 2 . Ej of the single 
junction SQUID depends on the magnetic field B, and 
it can be expressed as Ej — RqA/2Rn | cos(ttB / Bq) \, 
where Rn = 2.8 kfl is the normal state resistance of 
the SQUID, A « 200 AteV the superconducting gap 
of Al, and Bo — 16 G the measured period of the 
SQUID modulation. Here we assume that the two junc- 
tions in the SQUID are identical. Thus we calculate 
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FIG. 2: IVC of a single junction SQUID at two different values 
of magnetic field, corresponding to (A) Ej/Ec = 4.5 and (B) 
Ej/Ec < 0.2. At these values the biasing SQUID arrays had 
the same zero bias resistance _Ro ~ 10 MSI. The inset shows 
the two lowest energy bands of a single junction calculated 
for Ej/Ec = 2 (dashed line) and Ej/E c = 0.2 (solid line). 
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FIG. 3: IVC of a single junction in the overdamped qua- 
sicharge regime at different temperature (points) and the the- 
oretical prediction (solid lines). In the fit Vc = 30 fiV. 
From left to right the cryostat (fitting noise) temperatures are: 
T=50 mK (160 mK), 250 mK (260 mK), 300 mK (400 mK). 



Ej/Ec — 4.5 at the maximum (curve A). Experimen- 
tal uncertainties make determination much less accurate 
at the minimum of the SQUID modulation, where we 
estimate Ej/Ec < 0.2 (curve B). 

Both curves A and B in Fig. 2 show a Coulomb block- 
ade feature, followed by a back bending region due to 
the Bloch oscillations. The differences between these two 
curves are qualitatively understood from the Bloch band 
theory \Wi ■ The critical voltage, or maximum block- 
ade voltage, is determined by the shape of the lowest 
energy band, Vc — max[dE (g) / dq\. For Ej/E c > 1 
the lowest energy band becomes very flat (see inset of 
Fig. 2), which explains why curve A has a smaller block- 
ade voltage than curve B. Furthermore, the gap between 
the lowest and first excited energy band strongly depends 
on Ej/Ec (inset Fig. 2). The gap determines the max- 
imum current, or Zener break-down current Iz, above 
which Zener tunneling to higher bands leads to dissi- 
pation and increased voltage across the junction. For 
Ej/Eq = 4.5 the experiment shows Iz — 10 nA, whereas 
for Ej/Ec < 0.2, we measure Iz — 0.5 nA, in qualitative 
agreement with theory. 

This qualitative comparison can be made more quan- 
titative if we take in to account the effect of fluctuations 
due to the finite temperature electromagnetic environ- 
ment, which cause the measured blockade voltage to be 
smaller than the theoretical critical voltage Vc- Using 
duality arguments, Beloborodov et al. |ll| have calcu- 
lated the IVC of a junction in a high impedance envi- 
ronment with Ej/Ec > 1 in the presence of thermal 
fluctuations. In the theory the environment consists of 
a resistor (Z env — R) characterized by gaussian noise 
with temperature r no ; sc . Ref. assumes the noise to 
be white, i.e., to have classical correlations, which is the 



case as long as fc_BT no i sc > RqcVc/R. For our experi- 
ment, this condition is satisfied as long as T no i SO > 5 mK. 

In the experiment, the biasing arrays may generate 
more complicated fluctuations, but for the purpose of 
comparison we describe the arrays by one parameter, R, 
which represents the most basic way to characterize their 
fluctuations. For a bias current less than Iz, the IVC of 
the junction obtained in | ll| can be expressed in a form 
dual to that of H@ 



(V) = Vc3m 



Ii-ipehR/TviPeVc/n) 



I-i/3eI b R/7v(PeV C /n) 



(1) 



where [3 — l/A:sT no i se is the inverse of the noise temper- 
ature and I v {z) is the modified Bessel function of argu- 
ment z and complex order v. 

We have experimentally studied the influence of the 
cryostat temperature on the IVC of samples with the lay- 
out presented in Fig. IB. Here the single junction under 
test (area « 0.02 fim 2 ) is not tunable, and it is current bi- 
ased by two non-SQUID arrays consisting of 16 junctions 
with nominal area 0.01 /mi 2 . The IVC of these arrays 
(Fig. lb) shows a strong Coulomb blockade, providing 
a good current bias for the single junction. However the 
voltage probes are SQUID arrays, allowing the tunability 
of their impedance to an optimal value. We have found 
that the tunability of the voltage probes is essential for 
this type of experiment. 

In Fig. 3 the measured IVCs at fixed magnetic field 
(B = 262 G) and various temperatures are shown 
together with the theoretical curves calculated using 
Eq. (1). For these samples, the Josephson coupling of 
the single junction was large enough to realize the ex- 
act dual of the overdamped Josephson effect. In con- 
trast to previous works jl5j, we were able to measure 
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FIG. 4: Noise temperature as determined from fitting to the- 
ory (Fig. 3) plotted vs. the cryostat temperature. Inset: Mea- 
sured blockade voltage for three sets of data (symbols) and 
calculated blockade voltage for Vc = 30 fiV (solid line) if 

Jnoise — Tcryostat — T . 



IVCs showing a complete back bending to a zero voltage 
current, allowing for the first time a quantitative com- 
parison between the theory and experiment. Excellent 
agreement is achieved using one value of the parameter 
R = 150 kf2 for all curves, while adjusting the noise tem- 
perature T no i se for a best fit. For the theoretical curves 
the value Vc — 30 //V was used, which compares very well 
to the estimated critical voltage Vc = 28 ± 7 fxV. This 
value was calculated [l8j from the energies Ej = 270 fj,eV 
and Ec ~ 90 /ieV, determined from the junction param- 
eters as previously described (Rn = 2 kSl, C « 0.9 fF). 
The value R = 150 kQ, which represents the theoretical 
effective frequency independent impedance seen by the 
single junction, compares rather well with the measured 
i?o ~ 200 k£! of the voltage probes (SQUID arrays) at 
the particular magnetic field studied. These arrays ex- 
hibit a weaker Coulomb blockade than the non-SQUID 
arrays, hence the former determine the impedance which 
the junction sees. 



Figure 4 shows a comparison between T< 



cryostat 



and 



Tnoise j where the latter is determined by the fit to theory 
as shown in Fig. 3. In the region 175 mK to 250 mK, 
the two temperatures coincide, but the noise tempera- 



ture does not go below 175 mK. This saturation indicates 
the existence of residual noise with an effective temper- 
ature around 175 mK in the measurement system. As 
reported in [||, excess noise may be a consequence of 
inadequate filtering of the measurement leads. Another 
possible source of excess noise is shot noise in the cur- 
rent biasing arrays. These arrays have a strong Coulomb 
blockade, where charge transport is discrete and a shot 
noise is expected. We also find excess noise when the 
cryostat is above 250 mK. This temperature is close to 
the estimated odd-even free energy difference for the 
islands in the arrays, above which quasiparticlc tunneling 
in the arrays would lead to excess noise. 

Thus we see that noise reduces the measured blockade 
voltage below the theoretically predicted maximum crit- 
ical voltage Vc = 30 fiV. To further illustrate this point 
we plot the measured blockade voltage vs. Tcryostat in the 
inset of Fig. 4. The solid line is the theoretical predic- 
tion if T cr y OS tat = T no i so , independent of the parameter 
R. We see that a blockade voltage of the order of 10 /J.V 
would be expected at a temperature of 50 mK, if the ex- 
cess noise could be reduced. It is remarkable that this 
theory, which uses a minimal description of the fluctua- 
tions in the Josephson junction arrays, gives the correct 
order of magnitude for Vc(T), and very accurately re- 
produces the shape of the IVC. A more complex model, 
including frequency dependent impedance of the arrays, 
as well as additional sources of noise may give better cor- 
respondence, but at the price of many more parameters. 

In summary, we have experimentally studied the IVC 
of a single Josephson junction in a high impedance en- 
vironment (R > Rq) in the strong coupling regime 
(Ej > Ec), where the dual of the classical Josephson 
effect is realized. We show qualitative agreement with 
the Bloch band theory for the same single junction as 
the ratio Ej/Ec is tuned with the magnetic field. By 
taking in to account the finite temperature of the elec- 
tromagnetic environment, we can quantitatively explain 
the measured data with a minimal theory, where gaussian 
fluctuations are described by one parameter. 
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